Analysis of perturbed Boussinesq equation via novel integrating schemes

To analyze and study the behaviour of the shallow water waves, the perturbed Boussinesq equation has acquired fundamental importance. The principal objective of this paper is to manifest the exact traveling wave solution of the perturbed Boussinesq equation by two well known techniques named as, two variables (G′G,1G) expansion method and generalized projective Riccati equations method. A diverse array of soliton solutions, encompassing periodic, bright solitons, singular solitons and bright singular solitons are obtained by the applications of proposed techniques. The constraint conditions for newly constructed solutions are also specified. To enhance comprehension, the numerical illustrations of constructed solutions have been represented using surface plots, 2D plots and density plots. The results delineated in this paper transcend existing analysis, offering a novel, well-structured, and modern perspective. The solutions obtained not only enrich understanding of shallow water wave models but also exhibit efficacy in providing detailed descriptions of their dynamics.


Introduction
The world around us contains nonlinear phenomena and to describe these phenomena, nonlinear partial differential equations play a vital role.It is not possible to deny the importance of nonlinear partial differential equations (NLPDEs).NLPDEs have been widely used to explain and investigate the physical phenomena occurring in the world.Their applications have provided fruitful results in all fields of natural and social sciences including engineering and biological sciences [1].
In recent years, researchers have paid much attention to finding the traveling wave solutions of NLPDEs.Solitary waves and solitons are the type of traveling waves that were first discovered by J. Scott Russell in 1834.Solitons have some unique properties e.g.solitons come into existence when nonlinear and dispersive effects are canceled in a medium.Soliton acts like a single wave with one crust.It retains its shape while traveling at a constant speed.
Many researchers have worked hard for finding soliton solutions of many nonlinear evolution equations.The NLEEs are used to describe complex problems in various fields such as, plasma physics, mechanics of water waves, nonlinear optical fibers, control theory, describing oceanic and atmospheric influences, fluid mechanics, mathematical chemistry, biology, signal processing and many more [1,2].The methods of finding exact solutions of NLEEs using traveling wave transformations are a popular topic of research.Many useful methods have been employed for extracting solitary wave solutions, such as; improved tan �ðxÞ 2 À � expansion method [3], first integral method [4], modified auxiliary equation method [5], generalized exponential rational function (GERF) method [6], an inverse-G 0 G expansion method [7], the Lie classical method and unified method [8], the modified generalized Riccati equation mapping approach [9], the modified generalized exponential rational function method and the extended function method [10], new modified generalized Riccati equation mapping approach [11] and others.
The fundamental aim of this research article is to investigate perturbed Boussinesq equation (PBE) as it holds significant importance in fluid dynamics, particularly in the context of water waves.This model yields soliton solutions, shock waves, and singular solutions.These solutions are very useful in studying different wave behavior under diverse conditions.The perturbed quantities occurring in the model allow the reader to study more complex wave behaviors, including solitons and other solutions.The perturbed Boussinesq equation, which is used to describe the propagation of waves in shallow water is being investigated in this article.This equation incorporates various effects such as refraction, diffraction, shoaling, and weak nonlinearity in fluid dynamics.It is integrable equation and possesses soliton solutions.
Boussinesq equation in various formats has been analyzed by many scholars through different techniques such, as the extended hyperbolic function method [12], modified auxiliary equation techniques [13], G 0 G À � method [2] and approximate symmetry method [14].
The two variables G 0 G ; 1 G À � expansion method and generalized projective Riccati equations method are reliable and efficient approaches for obtaining new and novel solutions in the form of hyperbolic, trigonometric and rational functions.The extracted hyperbolic function solutions can further lead to bright soliton, dark soliton, Kink and anti-Kink solutions, singular solitons.The trigonometric solutions can be referred to as periodic solutions.All analytical methods have some limitations.The analytical techniques that are applied in this article also have some limitations.The proposed approaches extracted dark solitons, periodic waves and bright solitons.The approaches are strong and play an efficient role in finding solutions to variety of NLPDEs.The proposed techniques are widely employed in nonlinear dynamics and soliton theory, to create solitonic shapes.These techniques provide the full spectrum of soliton solutions.
The analytical methods that have been used in this article have significance importance in solving nonlinear partial differential equations.These methods have been used by many scholars to solve nonlinear equations, such as Konno-Oono equation [15], biological population model and KdV-Zakharov-Kznestsov equation [16] solved by two variables G 0 G ; 1 G À � expansion method.Whereas, generalized projective Riccati equations method is used to solve Klien-Fock-Gordon equation [17] and Lakshmanan Porsezian Daniel model [18].
This paper is divided into six sections: governing model is described in Section 2, Section 3 is about detailed explanation of above mentioned methods, Section 4 gives mathematical analysis of perturbed Boussinesq equation, Section 5 provides graphical illustrations and Section 6 concludes the whole paper.

Governing model (perturbed Boussinesq equation (PBE)) description
Boussinesq equation is the fundamental equation that is modeled for stability of waves by interaction of surface waves over shallow water waves.Boussinesq describes that this wave maintains its shape by balancing between precipitous effect of nonlinearity and smooth effect of dispersion.
The Boussinesq equation is given as where S(x, t) is a function of x and t representing surface elevation, α and ω are treated as constant parameters.
When ω > 0, then we end up with, a linearly stable, and the numerical computable equation Eq (1) known as Good Boussinesq (GB) equation.The solitary waves described by the GB equation solely occur for a finite range of velocities and can merge into one solitary wave [19].
If the sign associated to ωS xxxx is changed i.e ω < 0, Eq (1) becomes the well-known bad Boussinesq equation.It is used to describe a two-dimensional flow of shallow water over a flat bottom, assuming that the water waves have small amplitudes.Bad Boussinesq-type has its importance from both mathematical and physical points of view, but the research results on its initial boundary value problems are scarce.A part of the reason for the paucity of the results is due to the properties of the linear part of Eq (1) that are so bad that the traditional mathematical methods cease to be effective [20].
With the advancement in the Boussinesq equation, many new forms of Boussinesq equations are developed and are used frequently by researchers.In this respect, the perturbed Boussinesq equation (PBE) comes out, as where, ρ indicates stabilization term and ω represents coefficient of dissipation [21].PBE is a popular nonlinear evolution equation and it has been widely implemented in coastal, harbor, water and oceanic engineering.It helps in forecasting waves in coastal areas, breaking of waves, interaction of waves, shoreline circulation in intense and normal weather conditions [2].Moreover it is used to model tsunami waves, oscillations of tidal waves and also to model the characteristics of shallow water waves which occur at beaches, lakes and in rivers.PBE model has attracted the attention of many researchers.In [22] the perturbed Boussinesq equation is investigated using generalized Kudryashov method and sine-Gordon expansion method.A study is also conducted on new perturbed conformable Boussinesq-like equations to deduce soliton solutions [23].In [13], new approximate symmetry method is applied on (2+1)-dimensional perturbed Boussinesq equation and new soliton solutions have been derived.

Description of methods
This section contains the detailed description of two proposed analytical methods.

Method 1: Two variables
To have a full grasp of two variables G 0 G ; 1 G À � expansion method following preliminary points are important to keep in mind: Point 1.The linear ordinary differential equation of 2 nd order is considered, as Moreover, E and F obey the following ODEs where λ and μ are treated as constants.Point 2. When λ < 0, the general solution of Eq (3) is obtained, as where B 1 and B 2 indicate arbitrary constants.Consequently, F is taken as where When λ > 0, then the general solution of Eq (3) is obtained, as where F is considered as where the value of F has the form The NLEE is considered, as where R represents a polynomial in S(x, t) and its partial derivatives.Following are the important steps of two variable

Fist step. The following traveling wave transformation is used
Sðx; tÞ ¼ pðZÞ; where v represents a constant velocity.After applying this transformation, Eq (11) is converted to an ODE as where Q shows a polynomial of p(η) and its all derivatives with respect to η.

PLOS ONE
Second step.It is assumed that the Eq (13) has solution of the following form where c r (r = 0, 1, 2, . .., M) and d r (r = 1, 2, . .., M) are constants, which are determined later.Third step.In order to determine the integer M in Eq (14), the homogenous balancing is employed, through that the highest-order derivative and the nonlinear term occurring in the equation are compared.
Fourth step.By plugging Eq (14) into Eq (13) along with Eqs ( 4) and ( 6), the left-hand side of Eq ( 13) is transferred into a polynomial in E and F, in which the degree of F is not greater than 1.Setting each coefficients of gained polynomial to zero.The homogenous system of algebraic equations is obtained.Upon solving the system, the values of c r , d r , v, μ, B 1 , B 2 and λ for the case of λ < 0 are determined.
Fifth step.Following the same procedure as explained in Fourth step., inserting Eq (14) into Eq (13) besides with Eqs ( 4) and ( 8) for the case λ > 0 and for the case λ = 0 inserting Eq (14) into Eq (13) with Eqs ( 4) and ( 10)), the values of arbitrary constants c r , d r , v, μ, B 1 , B 2 and λ are extracted.After substituting the values of arbitrary constants in Eq (14), the exact solutions of Eq (13) have been obtained which are expressed by trigonometric functions, hyperbolic functions or rational functions for λ > 0, λ < 0 or λ = 0, respectively.

Method 2: Generalized projective Riccati equations method
First step.According to Method 2, the general solution of Eq (13) has the form where a 0 , a r and b r are constants which are determined later.The functions A(η) and B(η) satisfy the following ordinary differential equations A 0 ðZÞ ¼ �AðZÞBðZÞ; where valid for values of i = ±1 where H and δ are nonzero constants.When δ = 0 and H = 0 then solution of Eq (15) has the following form Second step.In order to determine the integer M in Eq (15), the homogenous balancing is employed, through that the highest-order derivative and the nonlinear term of the equation are compared.

Extraction of solutions for the proposed model
This section provides the exact solutions of perturbed Boussinesq equation by applying the two variable expansion G 0 G ; 1 G À � -method and the generalized projective Riccati equations method.The obtained solutions may be effective in providing detail description of shallow water waves models.In order to apply these methods, the following traveling wave transformation is considered.
In this paper, Eq (2) is solved for n = 1.The perturbed Boussinesq equation is converted into following ordinary differential equation by employing the wave transformation (25), as Upon integrating Eq (26) twice and taking the constants of integration to be zero, Eq (26) takes the form Application of both methods on Eq (27) is discussed in the following subsections.

Application of method 1 for PBE
This subsection is devoted for the application of two variable expansion G 0 G ; 1 G À � -method on PBE.Homogenous balancing principle gives M = 2. Substituting M = 2 in Eq (14), the following form of solution is obtained.
Putting Eq (28) into Eq (27) and utilizing Method 1, the following solution sets for two cases of λ have been derived.4.1.1Case 1: λ < 0. Hyperbolic solutions will be obtained in this case.The following solution sets are extracted by applying Method 1.
1 st solution set ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi 2 nd solution set 3 rd solution set The hyperbolic solutions extracted by 1 st solution set are obtained, as 3ðh À rÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi Particularly, by taking B 1 = 0 and B 2 = 1, Eq (29) yields the following solutions, as Again by setting B 1 = 1 and B 2 = 0, Eq (29) gives the solutions, as S 1 ðx; tÞ ¼ À 2lðh À rÞ j À 3ðh À rÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The hyperbolic solutions extracted by 2 nd solution set are obtained, as Taking B 1 = 0 and B 2 = 1, Eq (32) produces the following hyperbolic solutions, as Again taking B 2 = 0 and B 1 = 1, Eq (32) produces the following hyperbolic solutions, as The hyperbolic solutions extracted by 3 rd solution set are obtained, as Setting B 1 = 1 and B 2 = 0, Eq (35) gives the solutions, as Taking B 2 = 1 and B 1 = 0, Eq (35) gives the solutions, as 4.1.2Case 2: λ > 0. Trigonometric solutions will be obtained in this case.The following solution sets are extracted by applying Method 1.
ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi l 2 s 2 À m 2 q j ffi ffi ffi l p ; c 2 ¼ À 3ðh À rÞ j : 3 rd solution set The trigonometric solutions extracted by 1 st solution set are obtained, as S 4 ðx; tÞ ¼ À 3ðh À rÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi Taking B 1 = 1 and B 2 = 0, Eq (38) produces the following periodic solutions, as S 4 ðx; tÞ ¼ À 3ðh À rÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi Again taking B 2 = 1 and B 1 = 0, Eq (38) produces the following periodic solutions, as S 4 ðx; tÞ ¼ 3ðh À rÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The trigonometric solutions extracted by 2 nd solution set are obtained, as Setting B 1 = 1 and B 2 = 0, Eq (41) yields the following solutions, as Taking B 2 = 1 and B 1 = 0, Eq (41) gives the following solutions, as The trigonometric solutions extracted by 3 rd solution set are obtained, as Inserting B 1 = 1 and B 2 = 0, Eq (44) yields the following periodic solutions, as Again by inserting B 2 = 1 and B 1 = 0, Eq (44) yields the following periodic solutions, as Remark 1: The case when λ = 0 can be performed in a similar pattern.For the sake of simplicity this case is omitted in this paper.

Application of method 2 for PBE
This subsection is devoted for the application of generalized projective Riccati equations method on PBE.Homogenous balancing principle gives M = 2. Substituting M = 2 in Eq (15), the following form of solution is obtained.
Putting Eq.(4.23) into Eq.(4.3) and utilizing Method 2, the following solution sets for two cases of � have been retrieved.

When
O ˆSolution sets for this case are as follows: The extracted hyperbolic solutions corresponding to 1 st set are determined as follows: For i = −1, ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The extracted trigonometric solutions corresponding to 1 st set are given as follows: For i = −1, ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The extracted trigonometric solutions corresponding to 2 nd set are given as follows:   Remark 2: It is important to mention here that the similar nature of trigonometric solutions can be obtained by taking � = 1 and i = 1.Moreover, by taking � = 0, the rational solutions can be obtained.For the sake of simplicity, these cases have been omitted in this paper.

Graphical illustration
Graphical illustrations of few of the extracted solutions are provided in this section.The solutions that are obtained above indicate wave structures of different forms.These different wave structures can be proven useful in nonlinear evolution models.Therefore, to understand these varying wave structures, graphical illustrations are considered to be very effective.

Conclusion
The perturbed Boussinesq equation has been studied in this article using two distinct methods: the two variables G 0 G ; 1 G À � expansion method and the generalized projective Riccati equations method.These methodologies have demonstrated their efficacy not only in the context of the perturbed Boussinesq equation but also across various other nonlinear partial differential equations.Few limitations and restrictions are associated with every analytical technique.Techniques applied in this research too obey some limitations.By following the given restrictions on the proposed techniques we are successfully extract bright solitons, dark soliton, periodic solutions, singular solitons and combo solitons solution for the governing model.The applied approaches are very efficient in obtaining new soliton solutions for variety of NLPDEs.These methods retrieved diverse form of solutions, containing trigonometric, rational and hyperbolic expressions.On comparing our finding with [23][24][25], it has been found that the results obtained in this paper using suggested methods are new and have not been reported earlier in literature.To enhance the understanding of the physical properties of these solutions we have explained the existence of these bright solitons, singular solitons, bright singular solitons, and periodic soliton solutions, through graphical representations employing appropriately chosen arbitrary parameters.Moreover, density plots have been meticulously provided to facilitate a deeper comprehension of the obtained solutions.It is noteworthy that the results presented in this study are novel, offering fresh insights into the dynamics of shallow water,